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Abstract. Some of the rapidly oscillating (roAp) stars, 
have frequencies which are larger than the acoustic cut- 
off frequency determined from stellar models with atmo- 
spheres based on an Eddington or Hopf law. 

As the cut-ofF frequency depends on the T(t) relation 
in the atmosphere, we have computed models and adia- 
batic frequencies for pulsating Ap stars with T{t) laws 
based on Kurucz model atmospheres and on the Hopf's 
purely radiative relation. 

We compare the values of the cut-off frequency de- 
rived from expressions for the potential from Vorontsov 
& Zarkhov (1989), from Gough (1986), and from the ap- 
proximation of an isothermal atmosphere. These mod- 
els predict a different reflexion efflciency for waves and 
hence marginally different values for the cut-off frequency. 
The frequency-dependent treatment of radiative transfer 
as well as an improved calculation of the radiative pres- 
sure in Kurucz model atmospheres increase the theoret- 
ical acoustic cut-off frequency for roAp stars by about 
200/iHz, which is closer to the observations. 

Since evolutionary effects significantly influence the 
acoustic cut-off frequency we restrict the comparison of 
our computations with observations to those two 'patho- 
logical' roAp stars for which more reliable astrophysical 
parameters are available, HD 24712 and aCir, and com- 
ment briefly on a third one, HD 134214. For a Cir we find 
models with Kurucz atmospheres which have indeed a 
cut-off frequency beyond the largest observed frequency 
and which are well within the Toff ~ L/Lq error box. For 
HD 24712 only models which are hotter by about 100 K 
and less luminous by nearly 10% than what is actually 
the most probable value would have an acoustic cut-off 
frequency large enough. HD 134214 fits our models best. 



however, the error box for T^s - L/Lq is the largest of all 
three stars. 

One may thus speculate that the old controversy about 
a mismatch between observed largest frequencies and the- 
oretical cut-off frequencies of roAp star models is resolved. 
However, the model atmospheres have to be refined by 
investigating NLTE effects, among others, and the obser- 
vational errors for the astrophysical fundamental parame- 
ters have to be reduced further before a definite conclusion 
can be drawn. For the latter, asteroseismology can provide 
an important improvement by determining the frequency 
splitting 1^0 = (2 ^^{dr / c))~^ which is sensitive to the 
evolutionary status of pulsating stars. 

Key v^rords: Stars: atmospheres - chemically peculiar - 
osciUations - individual: HD 24712, HD 128898, HD 134214 
- variables: roAp 
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1. Introduction 

For 5 out of 28 known rapidly oscillating magnetic chem- 
ically peculiar (CP2, Preston 1974) stars, the so-called 
roAp stars, the largest observed frequency exceeds the the- 
oretical acoustic cut-off frequency, which is determined by 
the outermost stellar regions. 

Waves with frequencies larger than the cut-off fre- 
quency are not well reflected towards the resonant cavity 
and decay in the atmosphere with a decreasing amplitude. 
It has been argued by Shibahashi and Saio (1985) that 
the cut-off frequency is largely influenced by the T(r) re- 
lation which requires a careful modelling of these layers. 
Extensive studies of the external layers and of the cut-off 
frequency have been carried out for the Sun (e.g., Gough 
1986, Balmforth & Gough 1990). Investigations of other 
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types of stars have been neglected until the discovery of 
roAp stars. 

Frequently, atmospheres in stellar models are based on 
an Eddington or Hopf law (e.g. Mihalas 1978, called here- 
after standard models). In such a case the T(t) relation 
is given by 

= (3/4). T,Vr + g(T)), 

where T is the temperature. Teg the effective temperature, 
T the Rosscland optical depth and (/(r) = 2/3 for the 
Eddington law. These T(r) relations suffer from two major 
approximations: convection is not included and radiative 
transfer is considered to be frequency independent (i.e. 
the grey case is assumed). CP 2 stars with masses between 
1.6 to 2.2 Mq have two very thin convection zones below 
the surface. The outer zone extends into the atmosphere 
and it is therefore necessary to include convection in the 
model atmospheres. 

The ATLAS9-code developed by Kurucz (1991) takes 
into account the effect of frequency-dependent opacities 
from both continuum sources and from about 58 million 
lines, hence including blanketing effects. It also explicitly 
accounts for convection and radiative pressure. 

In this paper we describe our computations of stellar 
models where atmospheres have been implemented which 
were derived from the Hopf law as well as from the more 
consistent Kurucz model atmospheres. 

2. Stellar models 

In stellar structure calculations one separates the in- 
ner part, where the diffusion approximation for radiative 
transfer is valid, from an outer part, where the atmosphere 
is interpolated using a T(r, Teg, g) law, with the gravity 
g. When convection is included in a model atmosphere, 
the transition between the inner part of a stellar model 
and the atmosphere must be located at an optical depth 
T = Tb ^ 10, where the diffusion approximation is still 
valid (Morel et al. 1994). 

Effects from a magnetic field are neglected in our stel- 
lar models. To discuss the influence of convection and 
opacities on a depth dependence of the critical frequency 
we compare our results obtained with advanced model at- 
mospheres (based on ATLAS9) with a much simpler grey 
atmosphere (Hopf) which is a very liniitc;(l description of 
a stellar atmosphere, but which is frequently assumed for 
stellar models. 

2.1. Atmosphere 

In order to describe the atmosphere we used the LTE 
Kurucz ATLAS9 code without the "overshooting option" 
(Castelli 1996) to calculate an interpolation table for 
T{t,Tcs, g). Model atmospheres with solar composition 
were computed for logg = 4.2 and for T^g ranging from 
7400 to 10000 K, with steps of 100 K, and no additional 
contribution to line opacity by microturbulence has been 



Table 1. Fundamental astrophysical parameters for models of 
1.8 Mq at different evolutionary stages, with an atmosphere 
derived from Kurucz atlas9 model atmospheres. 



age log(Teff) log(L/Lo) R/Rq logg 
(10' yrs) 



40 


3.9209 


1.0162 


1.550 


4.312 


100 


3.9198 


1.0264 


1.574 


4.299 


225 


3.9171 


1.0426 


1.624 


4.272 


325 


3.9140 


1.0561 


1.673 


4.246 


400 


3.9113 


1.0645 


1.711 


4.227 


500 


3.9067 


1.0772 


1.773 


4.196 



assumed. The gravity varies in our 1.8 Mq model (Tab. 1) 
from log 5 = 4.313 at the ZAMS to log 5 = 4.196 at an 
age of 500 Myr. However, the approximation of a constant 
logfif {only for the atmosphere.^ in evolutionary calcula- 
tions is justified as the T(r, Ted, logg) law is rather in- 
sensitive to small variations of logg. Gravity effects on 
the cut-off frequency along evolutionary tracks have been 
studied by, e.g., Shibahashi (1991) and shall be discussed 
in subsection 3.4. 

Generally it is adequate to compute Kurucz models 
with 72 layers, but the numerical accuracy is insufficient 
when derivatives are important, such as the temperature 
gradient. This problem was also mentioned by D. Katz 
and C. van't Veer (priv. comm.) and we have therefore 
computed all our atmosphere models with 288 layers. 

For stellar models with a simpler atmosphere based on 
the Hopf law where convection is not treated it is sufficient 
to restrict the computations for the atmosphere down to 
Tb = 2. The values of qir) at each optical depth are in- 
terpolated from the values given in Table (3.2) of Mihalas 
(1978). 

The outer boundary of solar model atmospheres is usu- 
ally set to the bottom of the chromosphere, at Text ^ lO"**, 
where the minimum of temperature occurs for the Sun. 
Since it is not known whether CP2 stars have a corona 
or not, we have tentatively defined the outer boundary 
of our atmospheres at Text = 10~® and will justify this 
choice later. The model atmospheres, however, were al- 
ways computed up to r = 10"^-*'^^, as in Kurucz (1993) 
model grids. The density at the stellar model boundary 
was fixed to Pext = 1.895-lO^^^g • cm""^, which is the value 
for the Kurucz model atmosphere with T^e = 8080 K at 
the given optical depth. 

2.2. Internal structure 

We have computed representative models for CP2 stars 
of 1.8 Mq with the CESAM code (Morel 1993 and 1997). 
These models have about 1600 mesh points and include a 
ZAMS sequence. The EPF equation of state (Eggleton et 
al. 1973) was used, an initial hydrogen content X = 0.7, 
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and a heavy-element abundance Z = 0.02. The most re- 
cent OPAL95 opacities (Iglcsias & Rogers 1996) with a 
bi-rational sphne interpolation (Houdek & Rogl 1996) 
were incorporated. Convection is described by the stan- 
dard mixing-length theory with a mixing length A = aHp 
(where a = 1.4 is the mixing- length parameter and Hp is 
the pressure scale height). Because the two external con- 
vection zones are very thin, the general properties of the 
models are insensitive to the value of a (Gough & Novotny 
1993, Audard & Provost 1994). 

We shall call a "Hopf model" a full stellar model con- 
structed with an atmosphere derived from Hopf 's law, and 
a "Kurucz model" a model with a Kurucz model atmo- 
sphere implemented. 

The evolution of a 1.8 M© Kurucz model is summa- 
rized in Tab. 1, and the main characteristics of Hopf and 
Kurucz models with an age of 500 • 10*' years are given 
in Tab. 2. As expected, they have similar effective tem- 
peratures of about 8066 K, and their radii and therefore 
luminosities are also very similar. The positions of the con- 
vection zones, which move very little during the evolution, 
are also given. 



Table 2. Characteristics of our 1.8 M© models with an age 
of 500 • 10® yrs. The Hopf model has an atmosphere derived 
from the Hopf law, while the atmosphere of the Kurucz model 
is derived from Kurucz's model atmospheres. The first lines 
give the effective temperature, luminosity and radius in solar 
units, and the characteristic frequency Slg = (GM / R^Y^^ . The 
second lines indicate the positions of the two outer convection 
zones in units of the stellar radius R*. The deeper one is lo- 
cated between the radii Tci and rc2, the outermost convection 
zone lies between rc3 and rd- The positions of these limits are 
expressed in logarithm of the pressure in the third lines. 



model 


Teff- (K) 

rci/R* 
log Pi 


L/Lq 
rc2/R* 
log Pa 


R/Rq 
log Pa 


ng(10-'=rad.s-') 
logPi 


Hopf 


8067 


11.94 


1.772 


356.9 




0.9915 


0.9946 


0.9988 


1.000 




5.96 


5.53 


4.595 


3.97 


Kurucz 


8066 


11.95 


1.773 


356.6 




0.9914 


0.9945 


0.9988 


0.9999 




•J. 96 


5. ">;') 


1.57 


3.9:i 



3.1. Theoretical background 

Our oscillation code integrates the 4*'* order system of 

equations governing the stellar nonradial adiabatic oscil- 
lations and takes into account the gravitational poten- 
tial perturbations (Unno et al. 1989). These equations are 
obtained by perturbing the equations of momentum and 
continuity. They do not include the effects of a magnetic 
field. 

Because the light received from stars other than the 
Sun is integrated over the visible hemisphere, essentially 
only modes of low degree i can be observed. We have 
therefore restricted our calculations to to 3, and we 
consider high radial orders n, which justifies a Richard- 
son extrapolation between frequencies derived from mod- 
els with about 1600 and 800 points (Shibahashi & Osaki 
1981). 

To calculate the critical and the cut-off frequencies we 

make the following assumptions. For low-degree modes, 
the displacement is essentially vertical, so that the hori- 
zontal component can be neglected. Consequently, we con- 
sider only radial modes and, because we investigate modes 
of high radial orders, we adopt the Cowling approxima- 
tion, i.e. neglect the perturbation of the gravitational po- 
tential. 

The equation of motion for adiabatic oscillations can 
be written (Unno et al, 1989) as: 



$ = 0, 



(1) 



where $ = p^^^c^div^, and ^ is the fluid displacement, 
c = (FiP/ p)^/^ is the sound speed, w = 2tti' is the angular 
velocity (with the pulsation frequency p),V is the acoustic 
potential, and the derivative is with respect to the radius. 
Eq. 1 shows that acoustic waves are reflected towards the 
interior and are well trapped, if is smaller than V^, 
whereas if is larger, the mode propagates into the at- 
mosphere dissipating mechanical energy which decreases 
the mode amplitude. The critical frequency, Vc, is defined 
as the frequency above which modes propagate outwards. 
The cut-off frequency, um, is the maximum value of the 
critical frequency encountered in the outermost stellar lay- 
ers. 

According to Vorontsov & Zarkhov (1989), the poten- 
tial for radial modes can be written as: 



^ 2dr 



2 

c I - 

r 



9_ 

n1 



^dc2 

2c2 dr 



3. Critical and cut-off frequencies 

In this section we comment on the definitions of the critical 
(y^ and cut-off (I'm) frequencies, report on our calcula- 
tions of these quantities and compare them with observa- 
tions. 



c 



9 



J_dc^ 
2c2 dr 



(2) 



where N is the Brunt- Vaisala frequency. Another formu- 
lation is proposed by Gough (1986), assuming the charac- 
teristic length of the eigenfunctions (dln(5r/r)/dr)~^ be- 
ing short compared to r, so that the problem corresponds 
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to a plane-parallel layer in constant gravity. Under these 
conditions, the potential reduces to: 

V^=wl, (3) 

where uj^ = c^/{AH'^) (1 - 2diJp/dr) is the critical an- 
gular velocity and Hp = — (dlnp/dr)^^ is the density 
scale height. Vorontsov & Zarkhov (1989) use the acousti- 
cal depth as the dependent variable, while Gough (1986) 
uses the radius. The approximation of an isothermal at- 
mosphere is also often adopted for calculating v^, and 
Hp = Hp is then constant. In this case Eq. (2) and (3) 
reduce to: 



Rapid density variations in the medium, such as occur- 
ring at the transition between the bottom of the radiative 
zone and the top of the convective zone, lead to an abrupt 
rise in the potential which constitutes a reflecting bar- 
rier to waves. But to some extent, waves penetrate and 
may be even tunnelling outwards (sec e.g., Christensen- 
Dalsgaard and Frandscn 1983, Balmfortli & Gough 1990). 
The thickness of the potential barrier, which determines 
the reflexion efficiency, is a function of the outer boundary 
conditions and of the potential V. We shall see in the next 
section that the different expressions for V predict a dif- 
ferent thickness of this barrier and therefore a different ef- 
ficiency of the wave reflexion. In the solar case, a potential 
barrier occurs at r ~ 10~^, because the temperature rises 
above these layers and the potential behaves as T~^/^. 
This temperature minimum is related to the boundary 
between the solar photosphere and chromosphere. 

As there is no clear observational evidence for a corona 
around CP2 stars, the existence of a temperature min- 
imum and therefore of a maximum of the potential is 
not known either (Balmforth & Gough 1990). Note how- 
ever, that a temperature minimum could be caused also 
by other phenomena than an outside corona. The con- 
cept of the acoustic cut-off frequency as the largest criti- 
cal frequency, however, requires a decrease of the critical 
frequency after a local maximum is reached in the same 
atmospheric layer where a temperature minimum is ob- 
served. First evidence of a corona in Ap stars has perhaps 
been recently discovered by Simon & Landsman (1997). 

In the next section we shall speculate on the acoustic 
cut-off frequency, provided that a temperature minimum 
at an optical depth of either Text = 10~^ (as for the Sun) 
or 10~^ exists. The latter depth was chosen, because layers 
at r = 10^*' arc optically transparent for all wavelenghts 
(Castelli et al. 1997). Higher layers are even more sensitive 
to NLTE effects which are not included in our models, 
while deeper layers become opaque in various wavelength 
ranges. Thus 10^^ is a reasonable compromise. 

We denote , and vi^^ as the critical frequen- 

cies V/27r derived from the equations (2), (3) and (4), re- 
spectively, and similarly and as the cut-off 



frequencies, i.e the maximum values of the critical frequen- 
cies encountered in the upper atmosphere. 

3.2. Kurucz versus Hopf models 

Figure 1 gives the temperature profile as a function of 
the logarithm of the optical depth for stellar models of 
1.8 Mq with an atmosphere based on the Hopf law and 
on Kurucz models, as well as the specific Kurucz AT- 
LAS9 model atmosphere for Tea = 8080 K. We see that the 
Kurucz model reproduces very well the specific Kurucz 
model atmosphere of the same Toff. In the atmosphere, 
the frequency-dependent treatment of radiative transfer 
increases the local temperature at a given radius, except 
in the very outer layers where the temperature of the Ku- 
rucz model decreases below the limiting value of the Hopf 
law (Fig. lb). 

The depth dependence of the temperature has a direct 
consequence on the gradient, the ionization of H and He i, 
on Fi, and on c. It directly affects therefore the critical 
frequency (Eq. 1). Fig. 2 shows the temperature gradient 
V = d log T/ d log P and the adiabatic temperature gradi- 
ent Vad as a function of the logarithm of pressure, for the 
same models as in Fig. 1 (only the outermost convection 
zone is shown; the deeper one occurs at r/i? ~ 0.995, i.e. 
at logP ~ 5.5, see Tab. 2). The superadiabatic tempera- 
ture gradient exhibits two close peaks, which are related 
to the transitions between radiation and convection and 
back to radiation at the borders of the thin convection 
zone. For r > 0.03, the local temperature of the Kurucz 
model is larger than that of the Hopf model (Fig. 1) and 
the peaks of the temperature gradient of the Kurucz model 
are shifted therefore towards the surface compared to the 
Hopf model. Similarly, the ionization zones of H and He I 
and the peaks of the adiabatic exponent Fi are also closer 
to the surface (Fig. 3). 



Fig. 1. Temperature profile of models with an atmosphere 

derived from the Hopf law (dotted) and from Kurucz model 
atmospheres (dashed). Both have an effective temperature of 
about 8080 K. In the upper panel (a), the solid line corresponds 
to the specific Kurucz model atmosphere with Teff = 8080 K 
and log 3 = 4.2, and, as expected, coincide with the dashed 
line. The lower panel (b) shows the outermost layers. 



Fig. 2. Temperature gradient V = d log T/ d log P for the mod- 
els with an atmosphere derived from Hopf law (dotted) and 

from Kurucz model atmospheres (dashed). Both have an effec- 
tive temperature of about 8080 K. The adiabatic gradient Vad 
is also presented. The solid line corresponds to the specific Ku- 
rucz model atmosphere with T^s = 8080 K and log 5 = 4.2. 



N. Audard et al.: The acoustic cut-off frequency of Ap stars 



5 



At T < 0.03, (logP < 3.9), the temperature gradi- 
ent is radiative (Fig. 2) and the temperature in the out- 
ermost layers of the Kurucz model is smaller by about 
1220 K than that of the Hopf model (Fig. lb). On the 
other hand, the adiabatic gradient and Fi are larger for 
the Kurucz model (Fig. 3). While Fi reaches a maximum 
at logP ^ 3.7 for the Hopf model and then decreases 
monotonically, it slightly increases till logP 2 for the 
Kurucz model before decreasing. This behaviour has a di- 
rect influence on the acoustic cut-ofi^ frequency, as will be 
shown below. Fi decreases below 4/3 for the Hopf model 
and tends to 1. The Eddington law has the same result 
and furthermore a similar tendency is found (G. Houdek 
1997, private communication) for an atmosphere derived 
from the semi-empirical solar model 'C from Vernazza et 
al. (1981). The tendency of Fi to decrease below 4/3 re- 
veals the failure of T(t) relations based on the Hopf law, 
on Vernazza et al.'s model, or similar, to describe properly 
the upper layers of stars which are more massive than the 
Sun. 



Fig. 3. Adiabatic exponent Fi = ((?lnP/91np)ad for stellar 
models with an atmosphere derived from the Hopf law (dot- 
ted), and from Kurucz models (dashed). 

The potential barrier formed by V2 (Eq. 3) at the tran- 
sition between the radiative and the convective zones at 
logT between 1.5 and -1 is larger than for Vi (Eq. 2, Fig. 4). 
The consequences of an efficiently reflecting boundary 
(represented by V2) and a running wave boundary (ap- 
proximated by Vi ) have been investigated in the solar case 
by Gabriel (1992) and result in different mode amplitudes. 
In the upper layers the agreement between the three ex- 
pressions for the potential is better for the Kurucz model 
(Fig. 4b) than for the Hopf model (Fig. 4a). The discrep- 
ancies at r < 10~^ are probably due to NLTE and nonadi- 
abatic effects where the main assumptions for the formal 
expressions of the potentials break down. 



Table 3. Acoustic cut-off frequencies i^j^j , I'f;^-^ and i>^' in 
/kHz (sec Eqs. 2, 3 and 4), for the Hopf and Kurucz models. 
The boundary of the model atmospheres is assumed either at 
Text = 10~* or at 10~®. Tlic tlircc last lines indicate the differ- 
ence A between the respective Kurucz and Hopf models. For 
Text = 10~®, we also give the cut-off frequencies for models 
with an artificially suppressed convection in the outer layers 
(sec the next subsection). 



Model 




Text: 10-* 


10^" 


(no convection) 


Hopf 


(2) 
'^M 


2139 
2159 
2231 


2324 
2354 
2232 


(2061) 
(2345) 
(2224) 


Kurucz 




2453 
2495 
2463 


2528 
2743 
2465 


(2568) 
(2748) 
(2470) 




Ai-li 


314 
336 
232 


204 
389 
233 





Table 3 gives the values of the cut-off frequencies , 
ly^j^ and for the Hopf and Kurucz models, computed 
for the potentials Vi, V2, and V3, respectively, and with 
the assumption that the atmospheric boundary is located 
either at Text = 10"^ or at lO^*'. The differences of the 
frequencies computed for the Hopf models relative to the 
Kurucz models are substantial and amount to up to 10% 
of I'M- We must, however, keep in mind that the results 
derived from the Hopf law must be considered with much 
caution in optically thin regions, i.e. at t ^ 10^"^. 



Fig. 5. a: (top) Critical frequency Uc^^ (Eq. 2) for the Hopf 
(dotted) and Kurucz models (dashed); b: (center) same for 
i>c (Eq. 3); c: (bottom) same for (Eq.4). 



Fig. 4. a: (top) Logarithm of the potentials Vi, V2 and V3 
(Eq. 2, 3 and 4, full, dotted and dashed lines, respectively), in 
the upper layers of the Hopf model; b: (bottom) same for the 
Kurucz model. 

The behaviour of the critical frequencies , and 

for the Hopf and Kurucz models are compared in 
Fig. 5. In all cases, the potential in the upper layers is 
larger for the Kurucz model (at t < 0.01). Note that the 

maximum value of Vc^^ , the acoustic cut-off frequency t'^'* , 
does not occur at the smallest optical depth, but slightly 
deeper in the atmosphere. 



We have also compared the frequencies of acoustic 
modes of low and high radial order derived from the Hopf 
and Kurucz models. The latter have frequencies larger 
by about 1 /xHz compared to their Hopf counterparts. Al- 
though it represents less than 0.01% of the frequency, this 
difference is larger than the measurement error of ground- 
based observations and for the asteroseismic space mission 
COROT (Catala et al. 1995). 

3.3. Effects of convection 

As we have demonstrated in the previous subsection, Ku- 
rucz model atmospheres, taking convection and frequency- 
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dependent treatment of radiative transfer into account, 
considerably improve stellar oscillation models in the up- 
per layers. To specify which of these effects has the largest 
influence on the acoustic cut-off frequency, we investigate 
Kurucz models with a mixing-length parameter a = 0, 
hence turning off convection, and compare them with Hopf 
models. For the internal structure models we have imposed 
the gradient to be radiative throughout the envelope, and 
the convective core, of course, is not affected. Specific Ku- 
rucz ATLAs9-model atmospheres were also computed with 
a = 0. A comparison should therefore isolate the impor- 
tance of the radiative transfer relative to convection. 

In the models without convection ionization of H and 
He I occurs closer to the surface. As a result, the variations 
of Fi and of the temperature gradient also occur closer to 
the surface (Fig. 6). Since the outermost stellar layers are 
radiative the physical quantities are unchanged for mod- 
els without convection (Fig. 6) compared to models with 
convection (Fig. 2). Consequently, the cut-off frequencies 
have almost the same values (see Tab. 3 for the case of 
rext = 10-«'). 



Fig. 6. Temperature gradient of models with an atmosphere 
computed with the Hopf T(r) law (dotted), and with Kurucz 
models (dashed). The mixing- length parameter a was set to 
in Oder to suppress convection. The solid line represents 
the gradient of the Kurucz ATLAs9-model atmosphere with 
Teff =8080K. 

These results clearly indicate that the major factor in- 
fluencing the acoustic cut-off frequency is the inclusion of 
a frequency-dependent treatment of radiative transfer and 
a better calculation of the radiative pressure rather than 
the treatment of convection in the atmosphere. This result 
was indeed expected since the outermost layers, which are 
relevant for the cut-off frequency, are radiative. 

S.4- Effects of evolution 

The cut-off frequency also depends on the evolution, as 
was shown for example by Shibahashi (1991). From the 
ZAMS to the end of the main sequence and for models 
of 1.8 M0 this quantity decreases from 2855 to 1212 ^Hz 
for the Hopf models and from 3116 to 1328 /xHz for the 
Kurucz models. However, the relative difference between 
the cut-off frequencies from the Kurucz and Hopf mod- 
els remains almost unchanged along the main sequence 
and is about 8.5% (see Tab. 4). The cut-off frequency also 

1/2 

scales as the characteristic frequency fig = [GM/R^) 
(or as the frequency spacing, see Shibahashi 1991). We ob- 
tain similar quantitative results for all masses typical for 
roAp stars. When investigating a large parameter space in 
mass, age and metallicity, it appears to first order to be 
sufficient to compute models with the simple Hopf law and 



Table 4. Acoustic cut-off frequencies v'^ (Eq. 2) as a function 

of age for Hopf and Kurucz models with 1.8 Mq and Z = 0.02. 
The relative difference A is Kurucz-Hopf. 



age (Myr) 


Hopf 


Kurucz 


A in % 




(/xHz) 


(/.Hz) 




100 


2855 


3116 


8.37 


200 


2732 


2985 


8.48 


300 


2610 


2864 


8.86 


400 


2481 


2726 


8.99 


500 


2360 


2581 


8.56 


600 


2221 


2426 


8.45 


800 


1914 


2073 


7.66 


1000 


1510 


1650 


8.49 


1100 


1212 


1328 


8.76 



to approximate the effects of a better treatment of the at- 
mosphere on the acoustic cut-off frequency by increasing 
this frequency by about 7% to 9%. 

Since the uncertainty in the age determination can in- 
troduce a significant error in the computed acoustic cut-off 
frequency we will focus our discussion in the following sec- 
tions on those two roAp stars, HD 24712 and HD 128898, 
for which we have more reliable mass estimates due to the 
availability of hipparcos parallaxes. 

3.5. Comparison with observations 

For about 5 out of 28 known roAp stars, the largest pub- 
lished frequency (see Tab. 5 and references therein) ex- 
ceeds the expected theoretical acoustic cut-off frequency 
determined from standard stellar models. We do not con- 
sider here roAp stars for which the highest observed fre- 
quency probably is a harmonic of their nonlinear oscil- 
lation (HD 83368 (Kurtz et al. 1993), HD 101065 (Mar- 
tinez & Kurtz 1990), HD 137949 (Kurtz et al. 1991), and 
HD 161459 (Martinez et al. 1991)). 

As we have demonstrated in the previous sections, stel- 
lar models with a frequency dependent treatment of the 
radiative transfer in the atmosphere (Kurucz models) do 
have a higher cut-off frequency than models with a grey 
atmosphere (Hopf models) which used to be the baseline 
for most of the investigations in this field. This result is in 
agreement with speculations of Shibahashi & Saio (1985) 
and of Matthews et al. (1990, 1996) that a steeper than 
solar temperature gradient would increase the cut-off fre- 
quency and hence bring theoretical results closer to the 
observations. 

The photometric and spectroscopic properties of 
the roAp star HD 24712 with Tog = 7250 ± 150 K 
(Ryabchikova et al. 1997), log(L/L0) = 0.91 ± 0.04 
(based on tthipparcos = 0'.'02041 ± 0'.'00084, a bolomet- 
ric correction of -0.085 (Schmidt-Kaler 1982) and ne- 
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Table 5. List of roAp stars for which the largest published 
frequency is larger than the cut-off frequency determined from 
standard stellar models. 



HD 


Observed 


Ampl. 


Comment 




Vmax (/UHz) 


mmag 




6532 


2 402 


0.77 


Kurtz et al. (1996) 


24712 


2 807 


0.20 


Kurtz et al. (1989) 


128898 


2 566 


0.12 


Kurtz et al. (1994) 


134214 


2 950 


3.40 


Kurtz et al. (1991) 


203932 


2 838 


0.17 


Martinez et al. (1990) 



glecting interstellar extinction), can be reproduced with 
a model of 1.63 Mq, Z = 0.02 and an age of about 
900 Myr. At Text = lO"*^, the Hopf model gives a cut- 
off frequency = 2280 /iHz whereas the Kurucz model 
gives 2480 ^Hz. A stellar model with a Kurucz atmosphere 
hotter by about 100 K, less luminous by nearly 10%, and 
less evolved by 100 Myr would have a cut-off frequency in 
agreement with the largest observed frequency. However, 
such a model is compatible only with the lower left corner 
of the error box (see Fig. 7a). 

We have also computed Hopf and Kurucz models with 
appropriate age for HD 128898 (aCir). For the first 
time, a Kurucz model atmosphere was calculated with 
an opacity distribution function specific to the compo- 
sition of aCir (Piskunov & Kupka 1997). Stellar mod- 
els with 1.93 Afg, Z = 0.03 and an age of 400 Myr, fit 
the observed values (Kupka et al. 1996) of Tcff = (7900 ± 
200) K and \og{L/LQ) = 1.11 + 0.01/ - 0.02 (based on 
tthipparcos = 0'.'06097±0'.'00058, a bolometric correction 
of -0.12 (Schmidt-Kaler 1982) and neglecting interstellar 
extinction). The cut-off frequencies for the Hopf and 
Kurucz models are 2346 and 2600 ^Hz, respectively, at 
rext= 10~^. The acoustic cut-off frequency computed for 
the Kurucz model is compatible with the largest observed 
frequency (see Table 5) and we can therefore conclude, 
that no discrepancy may exist for this roAp star between 
theoretical and observed cut-off frequencies. 

However, the cut-off frequency depends on the model 
input parameters and one must therefore account for un- 
certainties inherent to observations and modelling. For 
example, there are problems with the photometric cali- 
bration of fundamental parameters of CP stars. Models 
of different mass and age can fit the same star in the 
H-R diagram within the error bars. For HD 24712, e.g., 
the cut-off frequency varies from 2725 /iHz (Kurucz model 
with 1.60 Mq and 800 Myr) which is only 80/iHz short of 
the largest observed frequency, to an even lower value of 
2294 /iHz (1.63 M0, 1100 Myr). For a Cir, a model with 
1.90 and 600 Myr gives = 2329 ^Hz which would 
be clearly smaller than the largest observed frequency. 

A similar situation exists for HD 134214 which, un- 
fortunately, has a considerably larger error box due to the 



HIPPARCOS parallax of only 10.92±0.89marcsec and a T^s 

which could be estimated only from photometric indices. 

Evolutionary tracks and lines of constant cut-off fre- 
quency for Kurucz models are plotted in Fig. 7a together 
with the observational error boxes. Models with Z = 0.02 
were investigated for HD 24712 and HD 134214, and mod- 
els with Z = 0.03 for HD 128898. The same physics as for 
our previous 1.8 M© models was used (see Sec. 2). Errors 
on the mass and age determination are about 0.02 Mq 
and 200 Myr for HD 24712, and 0.02 Mq and 100 Myr for 
HD 128898. 

Since the cut-off frequency of Kurucz models is larger 
than for Hopf models of same age, lines corresponding to 
the Hopf models are shifted to higher effective tempera- 
ture and lower luminosity. For Z = 0.02 the Hopf lines for 
= 2300,2500 and 2800 /xHz almost coincide with the 
Kurucz lines for z/^^ = 2500,2800 and 3000 /iHz, respec- 
tively, and for Z = 0.03 the Hopf lines for i^^^ — 2300 and 
2500 /LtHz coincide with the Kurucz lines for = 2500 
and 2600 /iHz, respectively. 

Asteroscismology is a powerful tool for determining the 
evolutionary status of stars via the frequency separation 

z/Q = (2 / (dr/c)) ^ ~ iyn,e - i^n-i,i 
Jo 

(see e.g. Shibahashi 1991, and Kurtz & Martinez 1993). 
If vq can be measured, a more reliable estimate of the 
evolutionary status can be derived than with our classical 
approach, because no bolometric correction (determined 
for chemically 'normal' stars) and interstellar extinction 
(with large local differences) have to be considered. Fig. 7b 
shows lines of constant frequency spacing for the same 
models as for Fig. 7a. The observed value vq = 68 /tHz 
for HD 24712 (Kurtz et al. 1989) is consistent with our 
classically determined error box and indicates an effec- 
tive temperature which should be larger by about 100 K 
than what was obtained spectroscopically. There is a se- 
rious problem for a Cir, because the observed value of pq 
is 50 /tHz (Kurtz et al. 1994) which cannot be reconciled 
with the spectroscopically determined effective tempera- 
ture and/or the luminosity. One has to stress, however, 
that the amplitudes for the overtone oscillations relative 
to the mode with the largest amplitude are very small and 
of the order of only a few 0.1 mmag ! 

We stress also that the frequency difference between 
models with atmospheres computed cither with the Hopf 
law or with Kurucz models is comparable to the changes 
which are obtained when introducing physical processes 
such as convective core overshooting (see, e.g., Audard et 
al. 1995). 

4. Conclusion 

We have shown that the major impact for modelling pulsa- 
tion frequencies close to the cut-off frequency comes from 
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Fig. 7. HR diagramme for stars with 1.58 Mq to 1.69 Mq for 
Z = 0.02 (age up to 1000 Myr), and with 1.90 and 1.93 Mq 
for Z = 0.03 (age up to 700 Myr) (dashed hnes). The roAp 
stars HD 24712, HD 134214 and HD 128898 are indicated by 
circles and error boxes. Pull lines are lines of constant cut-off 
frequency for the Kurucz models, for 2300, 2500, 2800 and 
3000 ^iHz for Z = 0.02, and for 2300, 2500 and 2600 ^Hz for 
Z = 0.03 (a), 

and lines of constant frequency spacing vq = (2 J^{dr/c))~^ 
for the same models, from 80 to 65 /xHz for Z = 0.02, and 
from 70 to 60 ^Hz for Z = 0.03. For HD 24712 the frequency 
splitting vo = 68 ytiHz (Kurtz et al. 1989), and 50 yuHz (Kurtz 
et al. 1994) for HD 128898 (b). 

the inclusion of a frequency-dependent treatment of ra- 
diative transfer and of blanketing effects, as well as from 
a better calculation of the radiative pressure in the model 
atmospheres, rather than the inclusion of convection. At 
the main sequence, Kurucz model atmospheres merged 
with stellar models derived from the CESAM code in- 
crease the cut-off frequency by about 8.5 % relative to the 
value derived from the Hopf T{t) relation. 

For two roAp stars with probably the best available 
mass and luminosity estimates, HD 24712 and aCir, we 
find models with Kurucz atmospheres and with parame- 
ters in agreement with the observational error box which 
have a theoretical cut-off frequency larger than the largest 
observed frequency and hence are in agreement with obser- 
vations. One may thus speculate that the old controversy 
about a mismatch between observed largest frequencies 
and theoretical cut-off frequencies of roAp star models is 
resolved. 

We have assumed solar abundances for our models ex- 
cept for a Cir, which, however, has an abundance pattern 
that does not significantly change the model atmosphere 
from one with solar abundances. HD 24712, on average, 
is even less peculiar than a Cir. Assuming that Z = 0.03 
reflects only a surface composition for a Cir and Z = 0.02 
would be the better choice for a stellar model, one can 
expect a smaller mass to fit a Cir and a smaller cut-off 
frequency as deduced from Fig. 7. For other roAp stars 
with a more peculiar abundance pattern, the deviation 
from model atmospheres with solar abundances or mod- 
els with a scaled heavy element abundance will be larger 
(Gelbmann 1997). Abundant rare-earth elements, through 
blanketing effects, could decrease the surface temperature 
and thus increase the cut-off frequency. This frequency 
might also be affected by a chemical composition gradi- 
ent (Vauclair & Dolez 1990). Finally, improving CP2 star 
pulsation models requires also to account for the magnetic 
field (Dziembowski & Goode 1996). 

The low surface temperature in the outer atmospheric 
regions provides a potential well which reflects the acous- 
tic waves inwards. The existence of a maximum value of 
this potential, which defines the cut-off frequency, requires 



an increase of the temperature in even higher atmospheric 

layers. While such an increase is observed in the solar chro- 
mosphere towards the corona, the first direct evidence of a 
chromosphere in Ap stars has only been recently shown by 
Simon & Landsman (1997). However, in the very outer at- 
mospheric layers nonadiabatic and NLTE effects are prob- 
ably relevant and should be considered in the future. Al- 
ready at T ^ 10"" NLTE effects might lead to an increase 
of temperature (Mihalas 1978) which would decrease the 
acoustic potential, thus defining the upper limit of the 
critical frequency, and hence the cut-off frequency. Un- 
fortunately, NLTE atmospheres with a treatment of line 
blanketing as sophisticated as in atlasQ, for at least the 
main elements contributing to the opacity in the upper 
atmospheric regions, are not available in the foreseeable 
future. 
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